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This  paper  introduces  a  systematic  approach  for  the  derivation  of 
models  that  describe  the  behavior  of  composites  with  a  periodic 
structure.  A  technique  for  analyzing  the  accuracy  of  the  models  is 
addressed.  Although  the  theory  has  been  developed  for  the  general, 
n-dimensional  setting,  the  main  ideas  in  this  paper  are  presented  for  the 
1-dimensional  case  only,  for  reasons  of  simplicity.  The  effectivity  of 
the  approach  is  illustrated  with  numerical  examples.  Also,  the  approach 
has  the  potential  of  being  used  to  adaptively  select  the  best 
homogenization  model. 
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1 .  INTRODUCTION 

The  mathematical  analysis  of  the  behavior  of  composites  that  have  a 
periodic  structure  involves  the  solution  of  a  partial  differential  equa¬ 
tion  that  has  highly  oscillatory,  periodic  coefficients.  Although  the 
structure  and  properties  such  as  existence,  uniqueness,  and  regularity  of 
solutions  of  this  equations  are  well  known,  it  is  impossible  to  compute 
the  solution  in  a  "standard"  way  due  to  the  complicated  "local"  behaviour 
of  the  solution.  Therefore,  various  "averaging"  techniques  that  separate 
the  "local"  «ind  "global"  scales  have  been  devised  to  develop  approximate 
solutions.  These  techniques  fit  under  the  general  category  of 
homogenization,  which  refers  to  models  of  "equivalent",  and  perhaps 
fictitious,  homogeneous  materials  that  are  derived  by  averaging  the 
properties  of  the  constituent  materials  of  the  composite. 

Homogenization  is  widely  used  in  various  contexts,  and  has  a  long 
history.  One  of  the  first  papers  ([1])  to  address  this  question  Is  more 
than  150  years  old.  Today,  there  is  a  large  literature  dealing  with  the 
mathematical  and  engineering  aspects  of  composite- type  materials.  For 
example,  [2]  and  the  survey  paper  [3]  address  the  mathematical  aspects  of 
the  problem,  whereas  [4],  [5],  snd  [6]  are  oriented  more  towards 
applications. 

Many  models  of  composites  are  derived  with  the  expectation  that  the 
model  will  have  a  variety  of  applications.  In  general,  however,  the 
homogenization  should  depend  on  the  particular  goals  of  the  analysis.  A 
typical  situation  is  discussed  in  [7]  and  [8]  where  it  is  shown  that  an 
analysis  of  shear  waves  and  longitudinal  waves  in  a  periodic  medium 
requires  different  homogenizations,  or  models.  The  relative  merits  of 
various  homogenizations  are  discussed  in  [9]. 
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In  this  paper,  we  will  discuss  Che  mathematical  concepts  and  the 
numerical  aspects  of  an  homogenization  that  is  appropriate  for  a 
particular  application  and  that  can  be  the  basis  for  an  adaptive  model. 
For  the  sake  of  presenting  the  ideas  clearly,  we  shall  restrict  ourselves 
to  a  one  dimensional  problem.  The  generalization  to  the  n-dimensional 
case  is  available  in  [10],  [11],  and  [12]. 


2.  THE  MODEL  PROBLEM  AND  A  HOMOGENIZATION  CLASS 

We  will  study  the  solution  of  the  ordinary  differential  equation 

(1)  L(uh)  =  (a(i)  (x))  +  b(|)  uh(x)  -  f(x),  x  €  (—,»), 

where  a  and  b  are  real,  2it-periodic  functions  satisfying 

0  <  Sq  <  a(y)  <  a,  and  0  <  bQ  <  b(y)  <  b, 

for  — it  <  y  <  it,  and  h  a  (small)  positive  number  corresponding  to  the 
period  of  the  medium. 

In  order  to  discuss  properties  of  u*1,  we  need  the  following  function 
spaces: 


{complex-valued  functions  u  defined  on 


/"  (  !  1^2  <x)l2  .,UI  dx  < 

—  j-0  dxJ 
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where  v  is  a  real  number.  It  is  shown  in  [10]  that  there  is  a  positive 
number  vg  (ag,a,bg,b)  such  that  for  any  0  <  v  <  vg  there  exists  a 
unique  solution  of  (1)  in  h!v  for  every  f  e  H®v. 

Next,  we  introduce  the  formalism  of  an  homogenization  class.  In 
sections  4-5  we  discuss  several  methods  of  selecting  the  various  elements 
in  the  homogenization  class.  We  mention  here  that  this  formalism  is 
sufficiently  general  to  include  all  of  the  homogenizations  that  are  used 
in  practice.  The  homogenization  of  this  class  consists  of  (complex¬ 
valued)  2ic-periodic  functions  { <k  :  k  «  1 . K}  and  the  triples  of 

constant  (complex)  coefficient  differential  operators  k  = 

1,...,K},  along  with  the  solution  Vk  of  the  equation 

(2)  AkVk  =*  \  f  for  k  -  1,...,K. 

The  function  W^,  defined  by 

(3)  "£<*>  -  X 

k=l 

is  then  used  as  an  approximation  of  u*1,  the  solution  of  (1).  The 

dependence  on  h  of  wjj  is,  for  the  most  part,  an  explicit  dependence 

on  h  by  some  of  (J^,  Ak,  1^,  and  C^,  and  not  just  on  the  argument 

~.  In  the  succeeding  section  we  will  often  refer  to  the  differential 
n 

operators  Ak,  1^,  and  by  their  respective  Fourier  symbols  ak(t), 

Pk(t),  and  yk  (t),  which,  by  our  assumptions,  are  polynomials  in  t 
and  may  depend  explicitly  on  h. 

In  general,  {Vk:  k  ■  1,...,K}  exhibits  the  global  behavior  of  the 
composite,  whereas  {<l>k:  k  ■  1,...,K}  exhibits  the  local  behavior  of  the 


composite 
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The  operators  {A^,B^,C^:  k  *  1,...,K}  could  be  easily  generalized  to 
pseudo-differential  operators,  although,  for  simplicity,  we  do  not 
consider  this  case. 

3.  THE  FOURIER  TRANSFORMATION 

It  is  shown  in  [10]  that  the  solution  u^1  of  (1)  can  be  written  in 
the  form 

(4)  uh(x)  =•  —/  f(t)4>(£,h,t)eitX  dt 

where  f  is  the  Fourier  transformation  of  f  and  y  -*■  <t>(y,h,t)  is  a 
2u-periodic  complex-valued  function  that  satisfies  the  differential 
equation 

(5)  -e“ihty  (a(y)  ^  (*(y ,h,t)eihty))  +  h2b(y)*(y,h,t)  =  h2. 

The  exact  meaning  of  the  integral  in  (4)  is  given  in  [10].  Properties  of 
<|>(y,h,t)  are  studied  in  detail  in  [11].  It  is  shown  in  [11]  that 
4>(y,h,t)  is  analytic  in  h  and  t,  by  which  we  mean  that  there  is  a 
complex  neighborhood  of  in  such  that  about  each  point  in  this 

neighborhood  of  R  ,  the  function  (h,t)  -►  $(*,h,t)  can  be  expanded  in  a 
power  series,  convergent  in  Hper,  in  which  each  coefficient  is  an 
element  of  Hper.  In  [13],  it  is  shown  that  <p<y,h,t)  -  0(|t-1|)  as  |tl 
♦  Other  representations  of  u*1,  developed  in  a  different  context,  can 
be  found  in  [2]  and  [14]. 

Now,  applying  the  Fourier  transform  to  (2)  and  (3)  yields 

®  ^  6  ( t) 

w£(x)  -  /  f(t)  l  ~T“r  Yk<t)<L(J)  eicx  dt. 

*  /2x  k-1  VtJ  k  k  h 


(6) 
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The  operators  k  »  1,...,K}  could  be  easily  generalized  to 

pseudo-differential  operators,  although,  for  simplicity,  we  do  not 
consider  this  case. 

3.  THE  FOURIER  TRANSFORMATION 

It  is  shown  in  [10]  that  the  solution  u^  of  (1)  can  be  written  in 
the  form 

(4)  uh(x)  *  -  /  f(tH(£,h,t)eitx  dt 

✓2*  —  h 

where  f  is  the  Fourier  transformation  of  f  and  y  -*■  <j>(y,h,t)  is  a 
2it-periodic ,  complex,  valued  function  that  satisfies  the  differential 
equation 

(5)  -e_lhty  ^  (a(y)  (<|>(y,h,t)elhty))  +  h2b(y)<Ky,h,t)  *  h2. 

The  exact  meaning  of  the  integral  in  (4)  is  given  in  [10].  Properties  of 

$(y,h,t)  are  studied  in  details  in  [11].  It  is  shown  in  [11]  that 

i|>(y,h,t)  is  analytic  in  h  and  t,  by  which  we  mean  that  there  is  a 

2  2 

complex  neighborhood  of  R  in  C  such  that  about  each  point  in  this 
2 

neighborhood  of  R  ,  the  function  (h,t)  *  4>(*,h,t)  can  be  expanded  in  a 
power  series,  convergent  in  Hper»  in  which  each  coefficient  is  an 
element  of  Hper*  la  [13],  it  is  shown  that  <J>(y,h,t)  =*  0(|t-i|)  as  |t| 
-*•  ®.  Other  representations  of  u*1,  developed  in  a  different  context,  can 
be  found  in  [2]  and  [14]. 

Now,  applying  the  Fourier  transform  to  (2)  and  (3)  yields 

(6)  w£(x)  -~J  f(t)  l  Ji—  Y  (tU  (J)  eitx  dt. 

K  /2u  —  k«l  akU'  k  h 
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In  section  2  we  noted  that  any  of  4^,  a^,  8^*  and  yk  may  depend  on 

h.  Upon  setting 


(7) 


K  3k(t) 

qK(y,h,c)  .  _Tk(tHk(y) 


it  follows  from  (4)  and  (6)  that 


uh(x)  -  wjj(x)  =•  ~  /  f(t)(Kph,t)  -  Qk(|,h,t))eitX  dt. 

Thus,  the  error  that  arises  from  approximating  u^1  by  wj}  depends  on  the 
difference 

(8)  <j>(y,h,t)  -  Qg(y,b,t)  . 

Consequently,  for  a  given  h,  we  want  to  find  2u-periodic  functions 
{ cl»k :  k  “  !»•••>£}  and  constant  coefficient  differential  operators  {A^, 
B^,  C^:  k  *  1,...,K}  so  that  (8)  is  small  in  the  sense  that 

(9)  /  /  1 4>(y,h,t)  -  QR(y,t)  |  2  I  g(t)  |  2  dt  dy 

-it  -® 

is  small,  where  either  g  =  f  or,  considering  a  class  of  functions  f, 
g(t)  -  (1+t2)  n,  for  example.  In  particular,  the  quality  of  a  given 
homogenization  and  the  range  of  its  application  are  determined  by  (9). 


4.  THE  FIRST  HOMOGENIZATION:  EXPANDING  IN  POWERS  OF  h 

As  we  pointed  out  in  section  3,  $(»,h,t)  is  an  analytic  function 

of  h  and  t.  Consequently,  we  can  write 


L 
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(10) 


4>(*,h,t) 


J>0(»ft)  +  4>  ^  ,  t  )h  +  (t7(*,t)h  +  •  < 


The  functions  {$^(»,t):  j  =  0,1,...}  can  be  determined  by  expanding  (5) 
in  powers  of  h  and  substituting  (10).  Solving  for  the  first  few 
<>j(*,t)  yields 


"<t>0(y,t) 

80(t)  2 

At  +B 

(11) 

< 

^(y.t) 

-  it  gQ(t)x1(y) 

^*2(y»c) 

=  (it)2g0(t)x2>2(y)  +  80(t)x2,0(y)  +  e 

”  ,(it) 

where 

1  r%  d*l 

A  =  2^  J  U(y)  +  a(y)  -  (y))dy 
-x  y 

(12) 

1  * 

B  =>  j  b(y)dy; 


P2(it)  is  a  fourth  degree,  even  polynomial  in  it  that  has  real 
coefficients;  and  Xj  >  *2  0*  *2  2  are  rea^~va^uec*  functions  in  Hper 

satisfying 

-  (a(y)  <y))  -  fy  <y>« 

d  0 

‘  d7  (a(y)  (y))  “  _b(y)’ 

■  It  (a(y)  (y))  "  “a(yKl  +  (y))  +  ^  (a(y)xj(y)). 
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on  (— it,x)  and 


it  it  it 

/  Xj(y)dy  =  /  x2>0(y)dy  =  /  x2>2(y)dv 


o. 


Substituting  (10)  and  (11)  into  (4),  we  can  suggest  the  following 
homogenization  in  terms  of  the  homogenization  class  that  we  introduced  in 
section  2,  because  the  variables  y  and  t  can  be  separated  in  each 
4>^(y,t)  and  because  <t>j(y,t)  is  a  rational  function  of  t. 


<I>J  =  1 . 

*2 

4^3  =  X2j0» 

*4  =  *2,2’ 

and 

cl  3  “  1 . 

The  Fourier 

symbols  of 

the  differential 

operators  are 

ac^t)  = 

At2 

+  B, 

?l<t) 

=  1, 

r^t) 

=  1, 

a2(t)  = 

At2 

+  B, 

32(t) 

“  1, 

Y2^c^ 

=  hit, 

II 

/-N 

u 

N-/ 

m 

At2 

+  B, 

63(t) 

=  1, 

Y3(t) 

=  h2, 

a4(t)  =» 

At2 

+  B, 

P4(t) 

=  1, 

Y4(t) 

=  h2(it)2 

a5(T)  - 

(At 

2  +  B)2 

>  ^5  ( tty 

=  1, 

Y5(t) 

=  h2. 

This  homogenization  is  accurate  when  f  is  smooth  and  h  is  small. 

Also,  only  the  Fourier  symbols  Yjc(t)  for  k.  >  2  depend  explicitly  on 

h. 

Setting  K  =  1  (i.e.,  using  clj,  <ij(t),  £3j(t),  and  y^(t)  only) 

yields  the  classical  homogenization  that  is  studied  e.g.  in  [2],  [15],  but 
is  derived  there  by  different  means.  The  asymptotic  rate  at  which  wj) 
approximates  u^1  as  h  -*•  0  is  studied  in  [2]  and  [15]  for  K  =  1,2. 


We  will  now  illustrate  some  of  these  ideas  with  an  example.  Let 

,  .  .  it 

1,  -  n  <  y  <  -  j 

(13)  a(y)  =  JlO,  -f  <  y  <  f 

^1,  <  y  <  rt 

and 

(14)  b(y)  =  1. 

Figures  la,b,c,  and  2a,b,c  contain  graphs  of  Re  <J>(y,h,t)  and 
Im  4>(y,h,t)  as  a  function  of  t  >  0  [<J>(y,h,t)  =  d>(y,h,— t) ,  4>(y,h,t) 

<)>(-yfhtt)  ]  for  various  values  of  y  and  h.  We  see  that  <fr(y,h,t) 
does  not  depend  on  y  very  much  when  h  is  small  and  when  t  is  not 
very  large. 


t  — 

Figure  la:  h  *  1. 


Re 
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t - - 

Figure  2c:  h  ■  1/10. 


For  h  -  1,  y,  figures  3a  and  3b  compare  Re  $(-. 8it,h,t)  with 

Re  QK(-.8ir,hft)  when  $q  (K  ■  1)  and  when  +  <^(»,t)h  + 

2 

$2(*>t)h  (K  *  3)  are  retained  in  (10)  (recall  that  Re  ^  ■  0).  If  h 
■  1,  it  is  clear  that  it  Is  best  to  use  just  one  term  in  (10).  However, 
if  h  ■  V2  ,  then  adding  the  third  term  ^('.t)!:  ,  will  improve  the 
approximation  if  |t|  is  not  too  large;  this  leads  to  increased  accuracy 
when  f  is  sufficiently  smooth.  A  similar  situation  occurs  when 
comparing  Im  $  with  Im  Q^;  see  figures  4a, b.  As  h  decreases  to  0 
we  see  two  effects.  First,  if  f  Is  sufficiently  smooth,  it  will  be 
advantageous  to  include  more  terms  in  (10).  Second,  a  given  truncation  of 
the  series  in  (10)  will  provide  a  good  approximation  of  #(»,h,t)  for 
larger  and  larger  values  of  |t|;  l.e.,  the  smoothness  restrictions  on 

f  can  be  relaxed  as  h  ♦  0. 


Re<X>,  Q  — ►  Re<&,  Q 
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Figure  3a:  h  *■  1. 


Figure  3b:  h  *  V2  . 
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Im  G>, 
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Figure  4a:  h  *  1. 


Figure  4b:  h  -  V2  . 
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5.  THE  SECOND  HOMOGENIZATION:  EXPANSION  IN  POWERS  OF  t 

In  section  4,  we  expanded  $(*,h,t)  In  powers  of  h.  Because 
$(«,h,t)  is  also  analytic  in  t  we  can  write 

(15)  $(*,h,t)  =*  $Q(*,h)  +  it  i>1(»,h)  +  (it^^C*  ,h)  +  ••• 

or 

06)  ^(. *h  t')~  ”  $g(*»h)  +  ic  4>1(*.1»)  +  (it)24>2(*  ,h)  +  •••  . 

Substituting  (15)  or  (16)  into  (4)  will  specify  the  functions  4^  and  the 
differential  operators  Ak,  B^,  and  Ck  that  should  be  used.  Note  that 
in  this  case,  the  functions  4>k  will  depend  on  h,  whereas  the 
differential  operators  can  be  chosen  independently  of  h. 

The  approximation  based  on  (16)  is  preferable  to  the  one  based  on  (15) 
because  4>(*,h,t)  ~  0(|t|  *)  as  1 1 1  -►<*>.  With  a  and  b  given  by  (13) 
and  (14)  in  section  4,  figures  5a, b  compare  Re  <t> (—  - 8x,h,t)  with 
Re  Q^(-.8x,h,t)  for  h  »  1  and  h  =  V2  ,  using  an  increasing  number  of 
terms  in  (16).  Figures  6a, b  compare  Im  <K-.8x,h,t)  and 
Im  QK(-.8it,h,t) . 

Comparing  figures  3a, b  with  5a, b  and  figures  4a, b  with  6a, b,  we  see 
that  the  quality  of  both  of  these  homogenizations  is  about  the  same, 
possibly  with  some  preference  for  the  homogenization  based  on  (16). 

In  [12],  the  homogenizations  that  are  described  in  sections  4  and  5 
here,  are  analyzed  in  detail  in  the  n-dimensional  case.  There,  estimates 
of  <fr(»,h,t)  -  QR(*,h,t)  in  various  norms;  e.g.,  H°er  and  Hper,  are 
given. 
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6.  OTHER  HOMOGENIZATIONS 

Obviously,  our  approach  gives  rise  to  many  other  homogenizations. 

Each  one  is  motivated,  more  or  less,  by  the  desire  to  minimize  (9).  Here, 
we  mention  some  obvious  choices. 

1)  Since  $(*,h,t)  is  a  2it-periodic  function,  we  write 


(17) 


<j>(y,h,t)  =  l  ck(h,t)elky. 

k“-® 


Mt) 

Approximating  ct(h,t)  by  the  rational  function  — — r-,  and  using  only  K 
K  ak(t) 

terms  in  (17)  yields  an  homogenization  that  is  in  the  homogenization  class 
described  in  section  2. 


2)  Instead  of  using  the  functions  e*"ky  as  a  basis,  construct  a 
basis  by  applying  the  Gram-Schmidt  orthonormalization  process  to  the 
functions  {  $(• .h.t^) :  k  -  1,...},  where  t^  is  chosen,  perhaps 
adaptively,  in  an  appropriate  manner.  Substituting  this  orthonormal  basis 
for  {eiky:  k  -  0,±l,...}  in  (17),  we  can  proceed  as  we  did  there. 


3)  For  "each"  y  f  (~rc,x),  approximate  $(y,h,t)  by  a  rational 

MO 


function 


a  (t)’ 
7 


and  solve 


Vy 


Byf .  Then 


uh(x) 


Wy(x) 


for 


—  (mod  2x) 
n 


and  x  ?  R. 


7 .  BOUNDED  DOMAIN 

Let  Q+  ■  {x|x  >  0}  and  Q_  -  {x|x  <  0},  and  let  the  differential 
operator  L  be  defined  as  in  (1).  Suppose  that  (1)  were  defined  on  Q+ 
with  the  boundary  condition  u(0)  ■  0.  Then  a  solution  can  be  developed 
in  the  following  manner.  Extend  the  function  f  into  Q_;  select  an 


("outside")  function  fg,  whose  support  is  small,  close  to  the  origin, 
and  contained  in  Q_;  and  solve  the  equation 

L(u^)  =  f  +  cfg  on  R* , 

where  c  is  chosen  so  that  uh(0)  a  0.  Note  the  choice  of  c  and  fg 
do  not  affect  the  value  of  L(u^)  on  0+  since  fg  =  0  there.  In  two 
or  more  dimensions  the  technique  is  similar. 

This  is  the  technique  that  is  used,  for  example,  in  the  theory  of 

elliptic  equations,  where  the  "outside"  functions  are  Dirac  functions  and 

the  error  on  the  boundary  is  minimized.  This  method  appears  to  be 
effective  for  solving  harmonic  or  biharmonic  equations,  especially  when 
the  position  of  the  Dirac  functions  are  comuputed  during  the  minimization 

procedure.  See  e.g.  [16]  and  [17].  For  problems  in  two  or  more  dimen¬ 

sions,  this  technique  needs  to  be  experimented  with  because  the  solution 
exhibits  a  boundary  layer  (see  [18]). 

8.  NUMERICAL  ASPECTS 

The  procedure  we  have  outlined  in  this  paper  for  deriving  approxima¬ 
tions  of  u*1  through  approximations  of  $(»,h,t),  has  certain  advatantages 
from  the  numerical  point  of  view.  A  (single)  finite  element  program  can 
be  used  to  solve  for  $(»,h,t)  as  defined  by  (5)  (also  in  two  dimensions), 
where  h  and  t  are  input  parameters.  Because  (5)  is  defined  on  one 
"cell"  only,  it  is  possible  to  solve  for  $(»,h,t)  with  a  high  degree  of 
accuracy  by  using  a  sufficient  number  of  elements.  The  derivatives  of 
$(«,h,t),  with  respect  to  h  and  t,  can  then  be  computed  by  numerical 
differentiation.  This  avoids  the  need  for  deriving  the  analytic  expan¬ 
sions  by  hand,  which  is  virtually  impossible  to  do  for  higher  orders.  The 
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graphs  In  this  paper  are  based  on  results  that  were  computed  by  this 
method. 

The  rational  functions  that  were  introduced  in  this  paper  can  be 
computed  in  a  standard  way.  The  functions  a^,  0^,  and  Y^  can  then  be 

derived  from  these  rational  functions,  and  used  to  construct  the  differen¬ 
tial  operators  A^,  B^,  and  of  section  2.  Solving  for  in 

equation  (2)  can  be  done  by  either  of  the  finite  element  or  finite  dif¬ 
ference  methods.  Furthermore,  it  would  still  be  possible  to  solve  (2)  if 

and  were  pseudo-differential  operators. 
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